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TECHNICAL MEMORANDUM NO. 534,: 

THEORY OF FLAPPING FLIGHT.* 
By Alexander Lapp is eh. 

The announcement of the 1925 Rhoen soaring- flight contest 
contains, among other things, the effort of a prize for flights 
in muscle-propelled aircraft. Obviously the aviator must 
first construct his flying machine and then try his luck with 
it. 

Before attempting to construct such a flying machine, the 
aviator will first try to post himself theoretically on the 
possible method of operating the flapping wings. It is not 
possible to include all the aspects of this difficult problem 
in the present article* W e will, however, give, on the basis 
of a simpler case, a graphic and mathematical method, which 
renders it possible to determine the power required, so far 
as it can be done on the basis of the wing dimensions. 

We will first consider the form of the flight path through 
the air* The simplest form is probably the curve of ordi- 
nary wave motion (Fig* 3, at the top)* This curve would be 
produced, in flight at uniform speed, by alternately moving 
the wings up and down (Fig* l). No account is taken of the 
fact that a longer stroke (if the wing of a bird is taken as 

a pattern) is made by the outer than by the inner portion of 
* From "Flugsport, 11 Juno 17, 1925, pp. 246-253. 
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a wing. The wing is here assumed to move as a whole, without 
the twisting of any portion with relation to the rest. The 
case is therefore as simple mechanically as can he conceived. 
On the other hand, the method of operation of a wing working 
like the wing of a bird can be derived from the results ob- 
tained in the simpler case. Moreover, this case assumes still 
greater importance when it represents the converse of dynamic 
soaring flight in air moving in the form of waves. 

After finding the flight curve, we must next determine 
the change in the angle of attack while passing through the 
different phases of the wave* Since the angle of attack, how- 
ever, is directly proportional to the lift, we will, for sim- 
plicity, determine the lift distribution along the flight 
path. We will refer the air forces themselves to the horizon- 
tal axis of the flight path and resolve them into a horizontal 
and a vertical component (Fig* 2). We thus obtain 

y^S^v* (1) 
x = c x S 1; V s (2) 

For the nondi mens tonal coefficients c y and. c x , we obtain the 
values 

e t , ~ c a + c w tancp . (3) 

c x " °a " fcan H ) ~ o w ( 4) 

Cy and c x are generally represented, as in the case of c a and-'Cw 
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by the 100-fold values Cy and e x . The angle q> is the momen- 
tary difference in direction 'between the horizontal and the 
tangent to the flight path and is found by the differentiation 
of the equation of the flight curve* The general formula is: 

y - m sin x ( 5) 

in the scale of the actual motion of the wing: 



m sin ( 2 ^ X A 
\ T o 

Herein: v = flight speed (m/sec/) , 

T G -~ time taken for a complete stroke ( sec . ) , 
with return to original position; 

x = distance flown (m), measured on the 
ho r i z ont a 1 1 inc. 

Then tan? ■= m cos ( 2 J ^\ (7) 

With the help of these formulas, we can calculate the mo- 
mentary values of x and y, or the corresponding values of 
c x and Cy, and plot them along the flight path as running 
functions of x. 

Since the angle cp is generally very small > it is not nec- 
essary to compute c y , because the member c w tan 9 is smaller 
than the possible error in e a . We therefore put Cy = c a (ap- 
proximately) and thus determine the lift distribution directly. 
c Y and c x are positive, when there is a thrust or negative 
drag. In order to determine the best kind of lift distribution 
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we will consider three cases: 

Case 1 (Fig. 3, a) r when the lift is constant. No change 
occurs in the angle of attack with relation to the flight 
curve. 

Case 2 (Fig. 3,b), when the lift fluctuates between its 
maximum value on the down- stroke and its minimum value on the 
up- stroke. The angle of attack with relation to the flight 
curve is greater during the down-stroke and smaller during 
the up- stroke* 

Case 5 (Fig. 3,c), when the lift fluctuates still more. 
The angle of attack with relation to the horizontal remains 
nearly constant (o^ = - 3.3°). 

Further data may be obtained from the diagrams. We learn 
from the diagrams that there is no forward wing thrust in 
case 1, while there is in case 2 and still more in case 3. 

The diagrams were drawn for a wing with the Gottingen 

profile No* 433 and an aspect ratio of 1 : 10, the maximum 

value of the angle q> ( in passing through the middle position) 

o 

being 9 lllax = ±8 and tan ~ 0*14 = m. 

The equation for the flight curve is accordingly 

0.140. sin ( g 2 1 loJL 

\ T Q V/ 2TF 

Any desired values may be given t- and v, so long as it 
is simply a matter of the computation or graphic representa- 
tion of the coefficients c x and Cy 
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For the computation of the required power and of the 
flight speed, we must determine the mean values of c x and c-y 
by the integration of the hatched areas. In order to do this, 
we represent the polar curves by the following equation: 

c w =' f-gP^ + c a 2 - p c a + q (8) • 

This form of equation for a parabola agrees best with the 
form of ordinary polar s because it takes account of the fact 
that the minimum resistance or drag occurs at small positive 
values of c a - When computed (by way of example) for profile 
433, it gives the following values: 

c w = + 0.0128^ c a s - 0.0093 o a + 0,0151 

The agreement fully suffices for our investigation and 
within the limits ^c a ~ 1.2 to c a - - 0*2, as shown by Fig. 4. 
The mean values obtained from the integration are then 



c- = I °a 



•'X 



mean 2 



Crr - C„ (9) 

j -mean a mean 

In the expression for Cr , Ac P denotes the fluctua- 
1 -"-mean °- 

tion of the different values of g & about its mean value. 
Hence, if (as in Fig. 3,c), the value of c a fluctuates be- 
tween c q = 1.2 and o flM+ = - 0 = 2, then the mean value of 
^max ^ffiin 

*= 1 .20 + (-0.20 ) = Qt BQ and ACa = ±o.70 



c a 18 



c 



a 
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We might, however, represent the distribution of c 8 by 
the following formula; 

c a v = Ac a C0S ( 2 - — f (10) 
co x a v t q y / ^niean 

The expression 



f ~ S — + c.n s - p c P + q 



"mean 



15 tt / 

i.e«.> it is the value of c w corresponding to the mean value 

Of 

Thus we had, by way of example, in case 3 (Fig. 3,c) the 
following as sump t ions J 

Ac a * ±Q.70 m - tan? r ,„ x = 0.140 ( -#- + = 0.0446 

Vb TT / 

(#" = To) c a m =0'50 o Wffl = 0.0216 

from which we then derive 

c = %^ [0.140 - 0.70 - 0.0446] = 0.0216 

•^mean 2 J 

c x „ = 0.0165 ( c x = 1.65) 

moan 

The maximum forward motion occurs for 

A C = 



VI? TT /' 



whereby c x becomes 



c = m 5 / - _ v 

x max 8 f s + k ^ ° w mean Ui; 



Vb tt / 
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This maximum is, however, of only conditional importance, 
"because &c a , even for small values of m, exceeds the burble 
point or the upper and lower interruptions of the polar curve. 

For our case, it would be 

AG ^ ~ 0^92 ' la3 

For the mean values of c„ =* 0.5. the limits would be 

Co = 2.23 and c~ . =? - 1.23. which far exceed the attain- 
cl max a min 5 

able limits of the polar. 

The wing thrust would then be 

c x = - 0.0216 - 0.0334 ( c* = 3.34) 

*max O.o563 A 

hence more than doubled. 

The ideal form of a flapping wing must therefore be sought 
in a wing whose polar curve covers a very wide range, as the 
result of suitable devices (adjustable profile, slots or rotor 
on the leading edge) « 

The resistance generated by the lift in the up and down 
motion of the wings and which must be overcome by the driving 
force is the total power consumption in the vertical thrust 
of the wings. 

The course of the work diagram for case 3 is shown in 
Fig. 5. 

By adding the work of the down- stroke and that of the up- 
stroke^ we obtain the hatched area in the form of an ellipse, 
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one of whose axes is Ac a . 

Consequently, the mean air resistance to be overcome toy 
the driving mechanism, is 

P m =0-78 Ac a S ^ v 2 (13) 

The vertical component of the flight path is 

S 0 = 2 m Z (13) 

Hence the force expressed in HP, is 

N 0.000415 m Ac a S v 3 (14) 

wheretoy ^ - -rg, i.e., for normal atmospheric conditions. 
Not/, however, 



v = 4 S - 



c„ 

■mean 



and hence 



g 



N = 0.021 m g V.-?,, , (15) 

a m 



Ac a 3 

If, in the example of case 3 (Fig. 3,e), we put g- 133 
we have 



3.18 



y f = 3.18 S = 12 m a and then 



S 

v = 4 rj^ft = 18 m/oec 
The required power is < 
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■N - 0*000415 x 0.140 x 0.70 X 12 x 5832 



or 



3.18 



N ~ 0.021 X 0.140 x 122 




N - 2.8 HP. 



0.5C? 
0,70 s 



Thereby it is tacitly assumed that the coefficient of the mean 
drag c w , increased by the coefficient of the structural 
drag c gJ is equr 1 to the forward thrust. The flight is there- 
fore uniform only when the following condition is satisfied* 



In Fig. 6 j I and c^ are plotted against each other 



for different values of c * This diagram shows that the re- 
quired power is proportional to the increase in the drag* 

The dashed line me plotted with reference to the horse- 
power of a propeller-driven airplane having the same flight 
speed* The efficiency, in this case, was assumed to be very 
good (up to 0.8) , since the efficiency of the wing-flapping 
mechanism had hitherto been disregarded. Moreover, energy is 
consumed only as the result of the resistance or drag due to 
inertia* This drag is, at most, very small for the right mu- 
tual adjustment of the flapping parts. 

In Figs. 7 and Q 3 the aspect ratio and the weight are 

plotted against the horsepower. Fig. 7 shows that the effect 

b s 

of the aspect ratio is very small above — = 10 and any 
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further gain would be more than offset by the necessary increase 
in weight* The weight increase plays a much more important 
part (the same as for a kite) • 

The horsepowers thus obtained are naturally minimum values,. 
On the other hand, we must bear in mind that it doubtless gives 
still more effective wing motions. 

We have also disregarded the fact that accelerations and 
retardations are produced (both horizontally and vertically) 
by the fluctuations in the engine power and the forward thrust 
which affect cither the shape of the flight curve or, if the 
latter is maintained, the flight path with respect to the 
aircraft* 

The accelerations naturally depend on the mass of the air- 
craft, as well as on the rapidity of the strokes. The greater 
the mass, the less rapid the strokes can be. 

If we take nature as our pattern, it appears that, with 
increasing weight, the time for a complete stroke in an air- 
craft like our example, may be much greater without any un- 
favorable result. Any mathematical computation would carry 
us too faro 

Another question appears to be of more practical impor- 
tance, namely, what would be the effect of other shapes of the 
flight path? We will discuss this question and the problem of 
flapping flight after the manner of birds in our next article. 

We will content ourselves today with the knowledge ob~ 



s 
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taincd from our simple example, which (notwithstanding the ap- 
parent technical difficulties of flapping flight, as compared 
with the advantages of propeller-driven flight) seems to just- 
ify serious attempts to solve this problem. 

Translation by Dwight M. Miner, 
lational Advisory Committee 
for Aeronautics. 
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Figs. 1,2 & 4 
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Fig. 3 
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Figs. 5 & 6 




Fig. 6 
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Figs. 7 & 8 
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Fig. 7 
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Fig. 8 



